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ABSTRACT 

We quantify the presence of active galactic nuclei (AGNs) in a mass-complete (A/* > 5 x 10^°Mq) 
sample of 123 star-forming and quiescent galaxies at 1.5 < z < 2.5, using X-ray data from the 4 
Ms Chandra Deep Field-South (CDF-S) survey. 41% ± 7% of the galaxies are detected directly in 
X-rays, 22% ± 5% with rest-frame 0.5-8 keV luminosities consistent with hosting luminous AGNs 
(io.5-8kcV > 3 X lO^^ergs"^). The latter fraction is similar for star-forming and quiescent galaxies, 
and does not depend on galaxy stellar mass, suggesting that perhaps luminous AGNs are triggered by 
external effects such as mergers. We detect significant mean X-ray signals in stacked images for both 
the individually non-detected star-forming and quiescent galaxies, with spectra consistent with star 
formation only and/or a low-luminosity AGN in both cases. Comparing star formation rates inferred 
from the 2-10 keV luminosities to those from rest-frame IR-I-UV emission, we find evidence for an 
X-ray excess indicative of low-luminosity AGNs. Among the quiescent galaxies, the excess suggests 
that as many as 70% — 100% of these contain low- or high-luminosity AGNs, while the corresponding 
fraction is lower among star- forming galaxies (43% — 65%). Our discovery of the ubiquity of AGNs 
in massive, quiescent z ^ 2 galaxies provides observational support for the importance of AGNs in 
impeding star formation during galaxy evolution. 

Subject headings: galaxies: active - galaxies: evolution - galaxies: high-redshift - galaxies: star- 
formation - infrared: galaxies - X-rays: galaxies 



1. INTRODUCTION 

Massive galaxies dX z ^ 2 a re likely progenitors of to - 
day's massive ellipticals fe.g.. ^ van Dokkum et al.|[2010[ ). 
Combining broadband UV-MIR photometry with high- 
resolution imaging has shown that a relation akin to 
the Hubble sequence is in place at z 2, dividing 
massive galaxies into two main categories: star- forming 
galaxies with disk-like or irregular structures and qui- 
escent galaxies with little or no star formation and ex- 
tremely co mpact morphologies (e.g., Cim atti et al.| [2008 
Toft et al.li200 9: WiUiams et al. 2010; W uvts et al 
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Szomoru et al.l 2012). The quiescent galaxies form a sig- 
nifica nt fraction (30%-50%) of all massive Z'y 2 galaxies 
(e^ lKriek et al.1 120061: 1 Williams et"all [20091: [ibft et alj 
120091 ) ■ but how they got to this stage remains an open 
question. 

Active galactic nuclei (AGNs) are one of the preferred 
hypothesized mechanisms for arresting star formation, 
as radiation, winds and jets from the central accreting 
black hole can remove the gas or heat it so that con- 
traction cannot take place (see iFabianI 120121 for a re- 
view of obser vational evidence) . Semi-anal ytical simula- 
tions by, e.g.. iDi Matteo et al.l |2005) and ICroton et al.l 
(|2006f ). show how reasonable descriptions of the AGN 
feedback on star formation in elliptical galaxies can re- 
sult in galaxies which resemble those observed in terms 
of, e.g., their mass function. This could be the main 
mechanism leading the massive galaxies at z ~ 2 from 
the blue star-forming cloud through t he "green valley" 
and o nto the red and dea d sequence (jSchawinski et al.l 
120071 : iHopkins et 311120081 ). An AGN would be capable 
of heating or expelling the gas residing in the galaxy as 
well as the gas that is believed to flow in via cold streams 



and otherwise prov ide the dominant source of star for- 
mation at z 2-3 (De kel et al1l2009f ). To help test this 
picture for massive galaxies at z ^ 2 we need better ob- 
servational constraints on the incidence of AGN and the 
connection to star formation in these galaxies. 

A general concern when estimating the star forma- 
tion rate (SFR) by fitting the broadband spectral en- 
ergy distribution (SED) of a galaxy with stellar popula- 
tion synthesis models, is the relative dominance of AGN 
versus star formation. The dominant contributor to the 
mid-infrared (MIR) light is from UV-light reprocessed by 
dust, but the UV light can be emitted by either young 
stars or an AGN. At z ~ 2 it is not clear, without the aid 
of high spatial and spectral resolution, what is causing 
the observed MIR emission from massive galaxies; is dust 
in the central regions being heated by AGN activity, is 
dust across a larger region of the galaxy being heated by 
star formation, or is a combination of the two scenarios 
taking place? 

X-ray emission, on the other hand, does not suffer from 
strong dust obscuration. Observing in X-ray can thus lift 
the apparent degeneracy in interpreting the origin of the 
MIR light, and thereby help to give the true AGN frac- 
tion of a sample of galaxies. A galaxy dominated by AGN 
activity can be distinguished from a galaxy dominated by 
star formation by having a harder X-ray spectrum if the 
AGN is obscured by dust, or by simply having a very 
high X-ray luminosity, io.s-skeV- The most heavily ob- 
scured "Compton-thick" AGNs, with column densities 
-/Vh > lO^"* cm~^, might be missed by X-ray selection, 
but can instead be identified via an excess of MIR emis- 
sion over that expected from purely star- form ing galaxies 
(e.g.. lDaddi et al.|[2007t iTreister et al.|[200l . 
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The shape of the MIR spectrum can also reveal AGNs, 
as the intense nuclear emission re-emitted by dust leads 
to a power-law spectrum in the MIR. For this purpose, 
color cuts have been devised using th e Spitzer da.ta. (e.g., 
IStern et al.l [20051: iDonTev et al.|[20Tl) . While this tech- 
nique has the capability of detecting even Compton-thick 
AGNs, otherwise missed in X-ray, it has a low efficiency 
for low- to moderate-luminosity AGNs and its rob ustness 
has yet to be verified at 2; > 2 ([Cardamone et al.l 12008). 

In a galaxy where star formation is dominant, the to- 
tal hard band X-ray luminosity, j^2-iokcVi pSuXi be used to 
estimate the SFR (e g..lGrimm et "all26oa^^ et al.l 
[200I iLehmer et al.ll2010D . An AGN will reveal itself if 
the X-ray SFR inferred this way is much larger than the 
SFR derived from other tracers such as Ha, UV, or MIR 
luminosity. Also, if star formation dominates the radia- 
tion output, the Lx ^> SFR conversion will give an up- 
per limit on the SFR as h as been obtained for sub-mm 
galaxies (jLaird "eraIl l2OT0|) and for massive, star-forming 



galaxies at 1.2 < z < 3.2 (jKurczvnski et al.|[2QT2l ). 

At redshifts around 2, high-ionization lines in the rest- 
frame UV can be used as AGN indicators, but X-ray 
observation s remain a, more efficient way of identify- 
ing AGNs JL amarei 113 120101 ). Several studies of mas- 
sive z ^ 2 galaxies have therefore been made with 
the aim of uncovering A GN fractions , using the Chan- 
dra X-ray observatory. iRubin et al.l ()2004 ) performed 
a study of 40 massive (M* = (1-5) x 10"Mq) red 
{Js — Kg > 2.3) galaxies at 2: > 2 by analyzing a 91 ks 
Chandra exposure. Roughly 5% of these were found to 
host an AGN with intrinsic i2-iokcV > 1.2 x lO^'^ergs"^. 
Assuming that the stacked X-ray signal from the re- 
maining X-ray undetected galaxies in X-ray comes from 
star formation alone, they derived a mean SFR broadly 
consistent with the typical mean S FRs estimated from 
SED fits. [Alexander et al.l ()20lH) analyzed the 4Ms 
Chandra observation s of 222 star-formi ng BzK galaxies 
(M, ~ 10i°-10"Mg. lDaddi et aT1l2007[ ) in the Chandra 
Deep Field-South (CDF-S). 10% (23) showed X-ray emis- 
sion consistent with AGN dominance, of which 5% (11) 
were found to contain heavily obscured AGNs, 4%(9) 
to have luminous, unobscured AGNs, and 3 out of 27 
low-luminosity systems showed excess variability over 
that expected from star formation processes, indicating 
that at least some low-luminosity (rest-frame L2-iokoV < 
XO^Sgj-gg-i^ systems may contain AGNs. 

The aim of this paper is to determine the AGN frac- 
tion in massive z ~ 2 galaxies, and reveal any differences 
between galaxies characterized as quiescent or star form- 
ing. We address the matter by analyzing the X-ray emis- 
sion from a mass-complete (M* > 5 x IO^'^Mq) sample 
of 1.5 < z < 2.5 galaxies residing in the CDF-S. The 
CDF-S, observed for 4 Ms, is currently the deepest X-ray 
view of the universe and so provides the best opportu- 
nity to study high-z galaxies across a relatively large area 
(464.5 arcmin^). 

Following a description of the sample selection (Sec- 
tions [2|) and of the X-ray data and analysis (Section [3]), 
we determine AGN fractions from the X-ray spectra and 
compare the X-ray inferred SFRs to UV-I-IR estimates 
(Section 14.11) . Results are presented and discussed in 
Section [5l and in Section [6l we sum up the conclusions. 
Throughout the paper we assume a flat cosmology with 



ftm — 0.3, f^A — 0.7 and h — 0.72, and magnitudes are 
in the AB system. 

2. GALAXY SAMPLE 

We select our gal axies from the FIREWORK^ catalog 
(iWuy ts et alll2008f ). which covers a field situated within 
the CDF-S and contains photometry in 17 bands from 
the U band to the MIR. For this study we extract a 
mass-complete sample of M^, > 5 x lO^'^M p, galaxies at 
1.5 < z < 2.5 in a w ay siin ilar to that of iFranx et al.l 
(|2008[ ) and lToft et all (1200^1. We will u se spectr oscopic 
redsh ifts when available (jVanzella et al.l [2008: Xuc et ahl 
120 lit) , and photometric redshifts from the FIREWORKS 
catalog otherwise. In order to maximize t he signal-to- 
noise (S/N) on results from X-ray stacking (jZheng et al.l 
I20TI) and ensure a relatively homogeneous PSF across 
the Chandra field employed, we consider only galaxies 
that lie within 6' of the average Chandra aimpoint. 

We ado pt galaxy st e llar n iasses from the SED fi tting 
results of iFranx et all pOOl and iToft et all pOOl ). In 
short, the SED fit s were carried out with models by 
iBruzual fc Charlotl (|2003l ). choosing the best fit from 
adopting three different star formation histories (a sin- 
gle stellar population with no dust, an exponentially de- 
clining star formation history of timescale 300 Myr with 
dust, and a constant star formation history with dust). 
In cases where the spectroscopic redshift differed by more 
than 0.1 from the original FIREWORKS photometric 
redshift, we redid the SED fits in FASlEl using an expo- 
nentially declining star formation history with a range of 
possible timescales from 10 Myr to ~ 1 Gyr. As a quality 
parameter of the SED modeling, we demand an upper 
limit on the reduced xt of 10 on the best-fit model. The 
SED fits provide SFR estimates, but we will be using 
SFRs derived from rest-frame UV-l-IR luminosities (see 
Section l4?2l) as these include obscured star formation and 
are subject to less assumptions. From the observed pho- 
tometry, rest-frame fluxes in U, T^. J band and at 2800 A 
have been derived using InterResiEl (|Tavlor et al.ll200l . 

We divide the resulting sample of 123 galaxies into 
quiescent and star-forming galaxies using the rest-frame 
U, V and J (falling roughly into the observed J, K and 
IRAC 4.5 /.tm bands at 2; ~ 2) colors. Dust-free but quies- 
cent galaxies differ from dust-obscured starburst galaxies 
in tha t they obey the following criteria bv Williams et al.l 
l|2009[ ): 

(1) 
(2) 
(3) 



U-V >1.3 
V -J <1.6 
{U -V)> 0.88 



x(y -J) + 0.49 



The fraction of quiescent galaxies identified within our 
sample using this method is 22% ± 5%(27/12 3). This is 
rather low compared to the 30-50% found bv lToft et al.l 
(|2009h for the same redshift and mass limit, but under 
the requirement that the sSFR (=SFR/M*) from SED 
fitting is less than 0.03 Gyr^^. If applying this crite- 
rion, we would arrive at a fraction of 42% ± 7%(51/123). 
A possible reason for the discrepancy between the two 
methods may be that we are using the UV J criterion in 

^ |http: // www. strw. leidenuniv.nl/f Ire works/l 
^ Ihttp : //astro .berkeley. edu/~ mar iska/FAST .html 
http : //www . strw. leidenuniv. Ill/- ent/InterRest 
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Figure 1. Stacked, background-subtracted and exposure- 
corrected 20" X 20" images in the three energy bands (columns) for 
the individually non-detected (top row) and detected (bottom row) 
quiescent and star-forming galaxies. Images of the non-detections 
have been smoothed with a Gaussian of width 2 pixels. 

the limits of the redshift range in which it has so far been 
established. As the redshift approaches 2, the quiescent 
population moves to bluer U — V colors, possibly crossing 
the boundaries of Equation ([3]). However, we prefer the 
UVJ selection technique in contrast to a cut in sSFR, 
because rest-frame colors are more robustly determined 
than star formation rates from SED fits. 

3. X-RAY ANALYSIS AND STACKING METHOD 

The raw X-ray data from the CDF-S survey consist 
of 54 individual observations taken between 1999 Octo- 
ber and 2010 Ju ly. We build our analysis on the work of 
IXue et ahl (|2011[ ) , who combined the observations and re- 
projected all images. They did so in observed full band 
(0.5-8 keV), soft band (0.5-2 keV) and hard band (2- 
8keV), and the resulting images and exposure maps are 
publicly available 

We extract source and background count ra tes in the X- 
ray m aps for all galaxies using the method of lCowie et al.l 
(|2012f ). Source counts are determined within a circu- 
lar aperture of fixed radii: 0".75 and 1".25 at off-axis 
angles oi 9 < 3' and 6* > 3' respectively. Background 
counts are estimated within an annulus 8"-22" from the 
so urce, excluding n earby X-ray sources from the catalog 
of IXue et al.l mm . 

Each galaxy is classified as "X-ray detected" if detected 
in at least one band at > Sa significance, thereby creating 
4 subsamples containing 8 quiescent and detected, 19 
quiescent and undetected, 43 star-forming and detected, 
and 53 star-forming and undetected galaxies. 

While the X-ray detected galaxies can be analyzed in- 



dividually, we stack the X-ray non-detected galaxies in 
order to constrain the typical X-ray flux from these. 
Stacked X-ray images of the 4 subsamples are shown 
in Figure [1] with all galaxies aligned to their Kg band 
galaxy center positions from FIREWORKS. Representa- 
tive count rates and associated errors for these stacks are 
ca lculated using the o ptimally weighted mean procedure 
of lCowie et al.l ( 2012|) and tabulated in Table [T] together 
with S/N values. 



Count Rate ztlcr (10" 


-«s-i) in 


Full Band 


Soft Band 


Hard Band 


Q 0.71 ±0.96 


2.01 ±0.41 


-1.67 ±0.86 


(0.7) 


(4.9) 


(-1.9) 


SF 2.51 ±0.55 


1.08 ±0.24 


1.29 ± 0.49 


(4.5) 


(4.5) 


(2.7) 



Table 1 

Stacking results for Quiescent (Q) and Star-forming (SF) galaxies 
not detected individually in X-Rays (in Parentheses the 
corresponding S/N) 



The reliability of our chosen method for X-ray stack- 
ing and source count extraction was tested with Monte 
Carlo (MC) simulations. With 500 MC simulations of 53 
(corresponding to the number of star-forming galaxies 
not detected in X-rays) randomly selected positions hav- 
ing no X-ray detections nearby, we obtain a histogram 
of stacked S/N values that is well fitted by a normal 
distribution with a center at —0.02 ± 0.8, that is, con- 
sistent with no detection. Similar results were obtained 
using a sample size of 19, the number of quiescent non- 
detections. 

We quantify the hardness of the X-ray spectra using 
the hardness ratio, HR== {H ^ S)/{H + S)E where H 
and iS* are the net counts in the hard and soft band, 
respectively. Assuming a power-law spectrum, we also 
derive a photon index, with the online mission count 
rate simulator WebPIMM Hjl using a Ga lactic H I column 
density of 8.8 x lO^^cm-^ ([stark et al.lll992f) and not in- 
cluding intrinsic absorption. Whenever the S/N in either 
soft or hard band is below 2, we use the corresponding 
2(7 upper limit on the count rate to calculate a limit on 
both HR and F, leading to a limit on the luminosity as 
well. When neither a hard- nor a soft-band flux is de- 
tected with an S/N above 2, a typical faint source value 
of F = 1.4 ( Xue et al. 2011;) is assumed, corresponding to 
HR^ —0.3 for an intrinsic powerlaw spectrum of F = 1.9 
(|Wang et al.ll200l . 

We derive the unabsorbed flux from the count rate by 
again using the WebPIMMS tool, now with the F tied 
to the observed value of either the individually detected 
galaxy or the stack in question. With this method, a 



count rate of ^ 10" 
nearly 10~^^ergcm~ 
of F = 0.8. Finally, 



counts s~ corresponds to a flux of 
s^^ in full band at a typical value 
the fiux is converted into luminos- 



^ Hardness ratio has the advantage, in comparison to F, of avoid- 
ing any assumptions regarding the shape of the spectrum. 

^ The photon index is defined as the exponent in 
the relation between photon flux density and energy: 



dAf/d_B [photons c m~^ s~^ keV~^] oc g~ 

* |http : //mm2 . astro . psu . edu/users/niel/cdfs/cdfs- Chandra .html ^ http : //heas arc . nasa ■ gov/Tools/w3plmnis ■ html | 
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-^0.5-8 kcV 


HR 


1,1 PGGifir'fi't'lrtn 


(crgs^i) 






> 3 X 1042 


< -0.2 


Unobscured AGNs (TVh < IO22 cm-2) 


> 3 X 10*2 


> -0.2 and < 0.8 


Moderately obscured AGNs (1022 < Nji < 102* cm-2) 


> 3 X 10*2 


> 0.8 


Compton-thick AGNs (Nu > 102* cm-2) 


< 3 X 10*2 


< -0.2 


Star-forming galaxy 


< 3 X 10*2 


> -0.2 


Low-luminosity obscured AGNs or star-forming galaxy 



Table 2 



Classification scheme used in this work, with limits from lSzokolv et al.l II2004I ) . FWang et al.l ||200'4I) . and ITreister et all II2009I) 



ity in t he desired rest -frame bands using XSPEC version 
12.0.7 dA rnau dl [T99l . 

4. INTERPRETATION OF X-RAY PROPERTIES 
4.1. Luminous AGN Identification 

In the top panel of Figure [2] are shown the rest- frame 
X-ray fuU band himinosities, io.s-skcVj not corrected for 
intrinsic absorption, versus observed hardness ratio, HR, 
for ah X-ray detected galaxies as well as stacked non- 
detections. A typical detection limit on io.5-8koV has 
been indicated with a dotted line in Figure [5J calcu- 
lated as two times the background noise averaged over 
all source positions. A few detected galaxies are found 
below this limit due to these residing at relatively low 
z and/or in regions of lower-than-average background. 
Galaxies with io.s-skcv > 3 x lO^^ergs"^ are selected 
as luminous AGN, since star -forming galaxies a re rarely 
found at these luminosities (jBauer etall[200l . In ad- 
dition, we adopt the HR criteria in Table [H in order to 
identify obscured and unobscured AGNs. 

About 53%(27/51) of the X-ray detected galaxies are 
identified as luminous AGNs, with 22 moderately to 
heavily obscured (—0.2 < HR < 0.8) and 5 unobscured 
(HR< —0.2) AGNs. The rest have X-ray emission con- 
sistent with cither low-luminosity AGNs or star forma- 
tion. We do not identify any Compton-thick AGNs di- 
rectly, but three galaxies have lower limits on their hard- 
ness ratios just below HR = 0.8, thus potentially consis- 
tent with Compton-thick emission (see Table[2]). In total, 
we detect a luminous AGN fraction of 22%±5%(27/123), 
a fraction which is 23% ± 5%(22/96) for the star-forming 
and 19% ± 9%(5/27) for the quiescent galaxies. 

Stacking the non-detections results in average X-ray 
source properties that exclude high luminosity AGNs. 
However, the inferred limits on their HR are consistent 
with a contribution from low-luminosity AGNs, the im- 
portance of which cannot be determined from this plot 
alone. 

4.2. X-ray Inferred SFR 

From the rest-frame hard band luminosity, ^2-10 kcV it 
is possible to derive estimates of the SFR, as the number 
of high-mass X-ray binaries ( H MXBs) is proporti o nal to 
the SFR (Ran aUi et a"Lll2003f ). iKurczvnski et al.l ()2012f ) 
made a comparison of different SFR indicators, by ap- 
plying three different SFR conversions to 510 star- 
forming BzK galaxies at 1.2 < z < 3.2, select ed as having 
. p2-ink ev < 10^'^ erg s~-^ . While relations bv iPersic et al.l 
(|2004f ) and iLehmer et al.l ()2010D overestimated the true 
SFR (from rest-frame UV and I R light) by a factor of 
~ 5, the relation bv iRanalli et al.i (,2003. ) provided a good 




-0.5 0.0 0.5 

Hordness ratio 



1.0 



1000 



> 
=) 
+ 



100 



10 



1 


"1 ^1 

/ 

/ 

/ 






/ • • 

' • 




—A 


1— • • • . 






1 ; 




• , 




1 1 1 





1 10 100 1000 10000 

X-roy (2-10keV) inferred SFR [Mg /yr] 



Figure 2. Top: Lo.5-8keV vs. HR for all galaxies. Red: quies- 
cent galaxies. Blue: star-forming galaxies. Squares: stacks of non- 
detected samples. Filled circles: luminous AGNs (selected with the 
X-ray criteria indicated by a shaded area). Open circles: galaxies 
dominated by low-luminosity AGNs or star formation processes 
alone (selected with X-ray criteria indicated by a hatched area). 
Dashed line: separating obscured (Afjj 1022 cm— 2) from un- 
obscured AGNs. Dotted line: typical detection limit. Errorbars 
display 2a errors and limits are indicated with arrows. Bottom: 
SFRuv+iR vs. SFR2-iokeV (see the text). Dashed line: equality. 
Typical errors on X-ray detected galaxies are shown in the lower 
right corner. Symbols are as above. 
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agreement at 1.5 < z < 3.2. But, as pointed out by 
iKurczvnski et al.l ()2012D . all relations might lead to an 
overestimation due to contamination by obscured AGNs 
in the sample of SF galaxies. As we do not know the 
exact amount of obscured AGN contamination in the 
stacks, we ha v e cho sen to use the following relation by 
iRanalli et afl (|200l as a conservative estimate of the 
SFR: 



SFR 



■2-10 kcV 



= 2.0 X 10-'^° L 



2-10 kcV 



(4) 



with SFR measured in Mq yr~^ and i2-iokGV in ergs~^. 
Another reason for n ot using the more recent relation by 
iLehmer et al.l ()201Cl[ ) is that this relation was constructed 
for galaxies with SFR > 9 Mq yr~^ only, whereas a large 
part of our sample has very low SFR s as inferred from 
SED fitting (< lM0yr~^). Following IKurczvnski et al.l 
([2012), we use the observed soft band flux to probe the 
rest-frame hard band luminosity. The uncertainty on 
the SFR is estimated from the error on the observed soft- 
band flux together with a sys tematic error in the r elation 
itself of 0.29 dex, as given bv lRanalh et~all (|200l . 

For comparison, the "true" SFR is inferred from rest- 
fr ame UV and IR l i ght, S FRuv-i-tr , following the method 
of iPapovich eTall (|2006l ): 



5. RESULTS AND DISCUSSION 

5.1. Luminous AGN Fraction 

In total, 22 X-ray detected galaxies have emis- 
sion consistent with containing a luminous (rest-frame 
-^o.5-8kcV > 3 X lO'^^ergs"^) and obscured (10^^ < A'^h < 
10^"* cm~^) AGN, and a further five have emission con- 
sistent with a luminous unobscured (A'h < 10^^ cm~^) 
AGN. This leads to a luminous AGN fraction of the full 
sample of 22%pr7i5%(27/123) and of the detected galax- 
ies only, 53% ± 13%(27/51). The AGN fraction among 
both quiescent and star-forming galaxies, according to 
their X-ray spectra, is measured to be around 20% as 
Table [3] shows, meaning that AGNs in massive z ^ 2 
galaxies, even quiescent ones, are common, as proposed 
by Krick et al. (2009) who studied the near-IR spectrum 
of one quiescent galaxy. 





Quiescent (27) 


Star-forming (96) 


Luminous AGNs 
Low-luminosity AGNs 


5 

2 (det) 
12-19 (non-det) 


22 
19 (det) 
0-21 (non-det) 


Luminous AGN fraction 
Total AGN fraction 


19% ± 9% 
70%-100% 


23% ± 5% 
43%-65% 


Table 3 



SFRi 



UV+IR 



= 1.8xlO-i"(LiR + 3.3L28oo)/iG 



(5) 



X-Ray derived AGN Numbers and Fractions for Quiescent and 
Star-forming Galaxies, Divided into Luminous AGNs and 
Detected and Non-detected Low-luminosity AGNs 



where Lir is the total infrared luminosity and L2800 is the 
monochromatic luminosity at rest-frame 2800 A. L2800 
comes from the rest-frame UV flux, /2800A (see Section 
121), and in this context, the errors on /2800A ^'"^ negligible. 
We derive Lir from the observed 24 flux, /24;im, using 
redshift-dependent multiplicative factors, a(z), from the 
work of .Wuvts et al.. (,2008. ') . Section 8.2 in that paper): 

LlR[Lo,8-1000/jm] — lO"'-^'' • /24/jm (6) 

The errorbars on Lir derive from the errors on /24,jm 
and we fu rther assume an uncer tainty of 0.3 dex in the 
relation of lPapovich et al.l (|2006D as found byHiil (|2003[ ) 
when comparing to Ha- and radio-derived SFRs. The X- 
ray undetected quiescent galaxies are not detected with 
S/N > 3 in 24 iim either (except in one case) leading to a 
mean flux of 1.2 ± 1.7/iJy, of which we adopt a 2a upper 
limit in the following. 

In the bottom panel of Figure [51 SFR2-10 kcV is com- 
pared to SFRuv-fiR with the dashed line indicating 
equality. It is no surprise that nearly all of the individual 
X-ray detections have very high SFR2-10 keV as compared 
to SFRuv-i-iR, as we are largely insensitive to individual 
purely star-forming galaxies, given the detection limits 
in the top panel of Figure [2] and the criteria in Table [2] 

For the star-forming stack, the X-ray inferred SFR of 
71 ± 51Af0yr~^ is consistent with the IR-I-UV inferred 
SFR of 46 ± 33M0yr~^, whereas the quiescent stack 
shows an SFR2-iokeV of 92±65Mo yr"! (89± 17Mq yr"! 
when bootstrapping this sample 200 times) , well exceed- 
ing SFRuv+iR < SMoyr"^ 



This luminous AG N fraction is high when compared 
to the 5% found by IRubin et all (|2004[ ) (see the intro- 
duction), but this is likely a consequence of their 4cr de- 
tection limit of 1.2 X lO^^ergs"^ in rest-frame 2-10 keV 
being about twice as high as our limit of 5.5 x lO**^ erg s~^ 
(at r — 1.4), as calculated from the average back- 
ground noise in t he observed soft ba nd. Adopting the 
detection limit of IRubin et al.l ()2004f ) and requiring an 
S/N of at least 4, we reduce our fraction of luminous 
AGN to 8% ± 3%(1 /123), consistent with t he resuhs of 
IRubin et al.l (|200l . [Alexander eFall (poTTh . using also 
the 4 Ms CDF-S data, found a much lower X-ray detec- 
tion fraction of 21%±3% as compared to ours (53%±7%), 
and a luminous AGN fraction of only 9% ± 2%(20/222). 
We believe that the discrepancies have several reasons, 
the main ones being: (1) our use of a mass-complete 
sample, whereas the B zK selection technique used by 
[Alexander et al.l ([201 1[) includes galaxies down to ~ 
10^*^ Mq for which the total AGN fraction, assuming a 
fixed Eddington ratio, is expected to be lower above our 
detection limitl^ (2) our updated source count extrac- 
tion and stacking method leading to higher S/N, and 
(3) the use of F instead of HR in the A GN identifica- 
tion conducted bv i Alexander et al.l ([201 ID . For compar- 
ison, iTanaka et al.l ([2012h recentlv discovered a group 
of quiescent galaxies at z = 1.6 with only one "main- 
sequence" star-forming galaxy. This group differed from 
local groups in having a remarkably high AGN fraction 

* For a fixed Eddington ratio, and assuming that galaxy bulge 
mass increases with total stellar mass, the AGN X-ray luminosity 
is expected to scale with galaxy stellar mass according to the A/th— 
Afbuigc relation ( Haring fc Jlix 2004) 
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of 38^20 J consistent with our result, and which they in- 
terpret as possible evidence for AGN activity quenching 
star formation. 

5.2. Importance of Low-luminosity AGNs 

As seen in Figure [2l the X-ray-identified luminous 
AGNs in general show an excess in SFR compared to that 
inferred from IR-I-UV emission. Among the galaxies clas- 
sified as being dominated by either low-luminosity AGN 
or star formation, about ~ 90%(21/24) have SFR2-iokoV 
more than la above SFRuv+iR- 

Surprisingly, the quiescent stack also has a much larger 
SFR2-iokcV than SFRir+uv- Even when removing the 
marginally undetected galaxies with 2 <S/N< 3, the re- 
sulting SFR2-iokcV = 62 ± 19Af0yr~^ is still more than 
3cr above SFRir+uv- This discrepancy is only further 
aggravated if instead assuming the SFR— relation of 
ILehmer et al.l (j2010,) . If caused by obscured star forma- 
tion, we would have expected an average 24 /zm flux of 
90 /iJy for the individual galaxies, in order to match the 
lower limit on SFR2-iokcV- This is far above the upper 
limit of 3.4 /iJy for the stack, suggesting that the X-ray 
flux of this stack is instead dominated by low-luminosity 
AGNs, but that their contribution to the 24 /ini flux re- 
mains undetectable in the Spitzer-MlPS data. 

We derive a lower limit on the low-luminosity AGN 
contribution for this stack of 19 objects by constructing 
a mock sample of the same size and with the same red- 
shifts, but containing N low-luminosity AGNs with lumi- 
nosities io.5-8keV — lO'*^ ergs^^ just below the detection 
limit and 19 — N galaxies with Lo.5-8keV = lO'^^ergs"^ 
(all of them with F = 1.4). From random realizations 
of this mock sample we find that at least N = 12 low- 
luminosity AGNs are required to match the observed 
rest-frame 2-10 keV luminosity of the stack. Similarly, 
we find that by removing at least 12 randomly selected 
galaxies it is possible to match the low SFR predicted 
by IR-I-UV emission, though only with a small probabil- 
ity 1% of 200 bootstrappings). In the following, we 
will adopt 63% (= 12/19) as a conservative lower limit on 
the low-luminosity AGN fraction among quiescent X-ray 
non-detections, but we note that the data are consistent 
with all the quiescent galaxies hosting low-luminosity 
AGNs if the luminosity of these is assumed to be only 
io.5-8kcV = 7 X lO^'^crgs-i. 

Interestingly, the quiescent stack is only detected in 
the soft band, cf. Table [TJ whereas one might expect a 
significant contribution to the hard band flux from a low- 
luminosity AGN population as the one proposed above. 
The lack of a hard band detection can be explained by 
the fact that the sensitivity of the CDF-S observations 
drops about a fa ctor of six from the soft to the hard band 
(IXue et al.ll2"0Tlll . meaning that the low-luminosity AGN 
population must be relatively unobscured (F > 1, con- 
sistent with the value of 1.4 assumed here), as it would 
otherwise have been detected in the hard band with an 
S/N> 2. 

The SFR2-10 koV of the star- forming stack is consistent 
with its SFRuv+iRj meaning that a strict lower limit to 
the low-luminosity AGN fraction here, is zero. However, 
performing the test above with the same model param- 
eters, a maximum of 40%(21/53) low-luminosity AGNs 
is possible, before the X-ray inferred SFR exceeds the 
upper limit on SFRuv-i-iR- 



It should be mentioned that the 24 /xm fiux, espe- 
cially at these high redshifts, is an uncertain estimator 
of the total rest- frame infrared luminosity, i.e., the entire 
dust peakj ^ used in the conversion to SFR. As shown by 
iBelll (|2003( ). one should ideally use the entire 8-1000 /xm 
range, e.g., by taking advantage of Herschel data, which 
can lead to systematic downward correction factors up to 

2.5 for galaxies with Lir w IO^^Lq (similar to the in- 
ferred IR luminosities of our sample galaxies detected in 
24 /im, showing a median of Lir — lO^^'^L©) as demon- 
strated by Elbaz et al. (2010). However, using the same 
conversion from f24i,.m to Ltr as the one implemented in 
this study, iWuvts et al.l (|2011[ ) showed that the resulting 
LiR for galaxies out to z ~ 3 are consistent with those 
derived from PACS photometry with a scatter of 0.25 
dex. Hence, we do not expect the inclusion of Herschel 
photometry in the present study to significantly impact 
any of our results, and we leave any such analysis for 
future work. 

Table[3]gives an overview of the derived AGN fractions, 
both at high and low X-ray luminosity. Adding the num- 
bers of luminous AGN, X-ray-detected low-luminosity 
AGNs as well as the estimated lower limit on the low- 
luminosity AGN fraction among non-detections, we ar- 
rive at a lower limit on the total AGN fraction of 



/agn > 



27 + 21 + 0.6 -19 + -53 
123 



0.48 



(7) 



for all massive galaxies at z '--^ 2. While for the star- 
forming galaxies this fraction lies in the range from 
43%-65%, it must be 70%, and potentially 100%, for 
the quiescent galaxies. Using the upper limits on these 
numbers, a tentative upper limit on the total AGN 
fraction is 0.72. 



5.3. Contribution from Hot Gas Halos 

We have so far considered star formation and 
AGN activity as causes of the X-ray emission ob- 
served, but a third possibility is an extended hot gas 
halo as seen around many nearby early-type galax- 
ies in the same mass-range as our sample galax- 
ies (jMulchaev fc Jeltemall2010D and predicted /observed 
around similar spirals dToft et al.''2002l: 'Rasm ussen et al.l 
120091: 1 Anderson fc Bregmanll2 011: Dai et al. 2012t). AGN 
X-ray emission is expected to come from a very small, 
R < 1 pc, accretion disk surrounding the centra l 
black hole of the host galaxy (|Lobanov fc Zens"u^l2007| ). 
whereas very extended star formation in the galaxy or a 
hot gas halo surrounding it would lead to more extended 
emission. We investigate the possibilities for these lat- 
ter cases by comparing radial surface brightness profiles 
of the 51 individually X-ray detected galaxies out to a 
radius of 8" in both the stacked observed image and a 
correspondingly stacked PSF image. 

The profiles are calculated in full band only, because 
of the high S/N here as compared to the other bands 
(cf. Figure [1]) . For each galaxy, we extract the back- 
ground subtracted source count per pixel within 10 con- 
centric rings of width 0".8 around the galaxy center po- 
sitions from the FIREWORKS catalog, and each profile 
is normalized to the mean count rate in all rings. The 
same procedure is applied to the corresponding PSF im- 
ages, extracted with the library and tools that come with 
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Detected galaxies 




Distance from center ["] 

Figure 3. Comparison of the radial surface brightness profiles in 
PSF images vs. observed images in full band for all detected galax- 
ies. Black: stacked observed image centered on K-hand positions 
from FIREWORKS for all detected galaxies (solid line) and if ex- 
cluding the three X-ray brightest sources (dashed line). Grey line: 
stacked PSF image. Red: stack of halo model images, with /3 = 0.5, 
Rc = 2.5 kpc (dashed) and 13 = 1, Rc = 1 kpc (dot-dashed). 

CIAO0 allowing for extraction at the exact source po- 
sitions on the detector. We verified the robustness of 
using PSF models from the CIAO calibration database 
for this purpose, by repeating the above procedure for 51 
kn own X-ray bright point-like sources from the catalog 
of iXue et al.l (|20lH l. and confirming that the resulting 
mean profile was fully consistent with that of the corre- 
sponding model PSFs. 

As can be observed in Figure |31 the combined radial 
profile of our X-ray detected galaxies in the full band is 
consistent with the PSF of point sources stacked at the 
same detector locations. A Kolmogorov-Smirnov (K-S) 
test yields a statistic of 0.2 and a probability of 99.96% 
that the two profiles are drawn from the same parent 
sample. 

Omitting the three sources with io.s-gkeV > 
10'*^ergs~-^ (cf. Figure [2]) leads to a more extended 
profile (dashed black line in Figure [3]), but the result 
still shows a high corresponding K-S probability of 70%. 
We also compare the profiles while recentering the im- 
ages on the center of the X-ray emission as derived using 
WAVDETECTEB instead of the i^-band center positions 
hsted in the FIREWORKS catalog, in order to test for 
the impact of any off-set between the X-ray and opti- 
cal centroids. However, the recentering only resulted in 
small variations within the errorbars on the original pro- 
file. The same conclusions apply to the subsample of X- 
ray detected star- forming galaxies only, whereas the S/N 
was too low to perform a similar study on the quiescent. 
X-ray detected sample alone. 

For comparison, we also simulated the stacked pro- 
file of extended hot gas halos, using for the sur- 
face brightnes s profile the /3- model first introduced by 
[Cavaliere fc Fusco-Femianol (J976). With default param- 
eters of /3 = 0.5 and core radius Rc = 2.5 kpc, both 
taken from the study oi z < 0.05 early-type galaxies by 

^ http: //cxc ■ cf a .harvard ■ 6du/clao4 . 3/aJielp/mk psf .html| 

1" WAVDETECT is part of the CIAO software. 



lO'Sullivan et al.l (|2003D . the stacked profile of the halos, 
convolved with the corresponding PSFs, is well outside 
the errorbars of the measured profile. Only if assuming 
an extremely compact case of /3 = 1 and Rc — I kpc, does 
the halo emission become sufficiently compact to mimic 
that of the PSF profile (see dash-dotted line in Figure|3]), 
but the halo model then overpredicts the observed sur- 
face brightness on scales smaller than 1". In conclusion, 
a hot gas halo alone, described by a /3-model, cannot ex- 
plain the emission, unless one chooses model parameters 
that render the profile indistinguishable from a point-like 
source. 

5.4. Quenching of Star Formation by AGNs? 

It remains debated whether the presence of AGN 
is connected with internal galaxy properties associated 
with secular e volution or, to a h igher degree, with exter- 
nal processes (|Darg et al.ll2010| ). For example, the cos- 
mic star formation rate density and the number density 
of AGN share a common growth history with a peak in 
activity around z ~ 2-3 (jTreister fc Urrvir2012l : iDunlopI 
1201 Ih . hinting at a co-evolution betwe en SFR and super - 
massive black hole (SMBH) accretion (|Schawinskill2011| ). 
In this section, we will investigate the correlation, if any, 
between the presence of AGNs and internal / external pro- 
cesses, treating luminous and low-luminosity AGNs sep- 
arately. As typical internal properties governing secular 
evolution we will focus on stellar mass, M*, and SFR 
(Peng et al. 2010), while major mergers are taken as a 
likely case of external processes and as a phenomenon 
often associated with luminous AGNs (jTreister et al.l 
[2011 . 

Starting with our X-ray identified luminous AGNs, we 
already found that the fraction of these does not corre- 
late with star formation in our sample (cf. Table [S]) . In 
addition, the distribution of luminous AGNs and that 
of the rest with regards to their SFRir+uv are similar, 
with a K-S test yielding a statistic of 0. 2 1 wit h a proba- 
bility of 59%. Similarly, iHarrison et al.l ()2012| ) found no 
correlation between AGN luminosity and quenching of 
star formation, albeit at a higher luminosity (L2-8kcV > 
10'*4ergs~^) than probed here. Dividing the sample ac- 
cording to M* instead, by constructing two bins around 
the median mass of 1.1 x 10^^ Mq, we arrive at similar lu- 
minous AGN fractions above and below this mass limit, 
namely 21% ± 7%(13/61) and 23% ± 7%(14/62), respec- 
tively. We thus find no clear evidence for the luminous 
AGN fraction of our sample to correlate with internal 
properties, suggesting that an external factor is of larger 
importance. 

An alternative is that luminous AGNs are primarily 
trigge red by non-secular p rocesses such as major merg- 
ers. iTreister et al.l ()2012[ ) found that major mergers 
are the only processes capable of triggering luminous 
{Lho\ > lO^^grgs-i) AGN at < z < 3. Our luminous 
AGN fraction is consistent with the major merger frac- 
tion of m assive (M^, > lO^- '^Min) galaxies at 1.7 < z < 3 
found by iMan et al.l (|20l2l ) to be 15% ± 8% (compared 
to a luminous AGN fraction of 11 ± 3%(13/123) when 
excluding galaxies with Af, < IQ^^Mq). This is consis- 
tent with the idea that our luminous AGNs are triggered 
by major mergers, but a more direct test of this scenario 
would be to search for major mergers using the imaging 
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data available in CDF-S. Indeed, iNewman et alj (|2012j ) 
used iJSTCANDELS imaging in the GOODS-South field 
residing within CDF-S (along with the UKIRT Ultra 
Deep Survey) to arrive at roughly equal pair fractions 
when comparing massive (M* > lO^^'^M©) quiescent to 
massive star- forming galaxies at 0.4 < z < 2. Again, 
this is in agreement with our result that the luminous 
AGN fraction does not vary between quiescent and star- 
forming galaxies. A detailed morphological study of our 
X-ray identified luminous AGN using the HST/WFC3 
data to search for evidence of merging is an interesting 
extension to the work presented here, which we leave for 
a future study. 

Turning toward our low-luminosity AGNs, we do see 
X-ray evidence for an enhanced population among the 
quiescent galaxies when compared to their star-forming 
equivalents. The mean masses of our non-detected sam- 
ples of quiescent and star-forming galaxies are similar 
(14.6 ± 1.6 and 12.7 ± 1.0 x 10^° Mq respectively), sug- 
gesting th at the relevant fac t or he re is SFR. At 0.01 < 
z < 0.07, ISchawinski et all (|2009[) found that massive 
(M* > IO^^Mq) host galaxies of low-luminosity AGNs 
all lie in the green valley, that is, at some intermedi- 
ate state after star formation quenching has taken place 
and before the SED is truly dominated by old stellar 
populations. The observation that these host galaxies 
had been quiescent for ^ 100 Myr, ruled out the pos- 
sibility that one short-lived, luminous AGN suppressed 
star formation and, at the same time, made it unlikely 
that the same AGN quenching star formation was still 
active, given current AGN lif etime estimates of ~ 10''- 
10^ yr (|Di Matteo et al.ll2005D . Rather, the authors fa- 
vored a scenario in which a low-luminosity AGN already 
shut down star formation, followed by a rise in lumi- 
nosity, making the AGN detectable. At z ~ 2, SED 
fits of massive (M* > IO^^Mq) quiescent galaxies show 
that the quenching typic ally took place ~ 1 Gyr before 
the time of observation ()Toft et al.l 120121 : Krogager, J.- 
K. et al., in preparation), demanding an even longer de- 
lay or an e pisodic AGN activ ity as frequently applied 
in models (|Croton et al.l l2006f ). Episodic AGN activ- 
ity could explain why we see evidence for a higher low- 
luminosity AGN fraction among quiescent as compared 
to star-forming galaxies, but it would also require the 
low-luminosity AGN phase in quiescent galaxies to last 
at least as long as the dormant phase. Future modeling 
and observations will show whether this is in fact possi- 
ble. 

We conclude that our data are consistent with a sce- 
nario in which luminous AGNs in massive galaxies at 
z ~ 2 are connected with major mergers or other non- 
secular processes, while the presence of low-luminosity 
AGNs in the majority of quiescent galaxies suggests that 
these AGNs present an important mechanism for quench- 
ing star formation and keeping it at a low level. Ulti- 
mately what happens at z ^ 2 has to agree with the 
subsequent evolution that changes size and morphology 
of quiescent galaxies (Barro et al. 2012) 

6. SUMMARY 

Our main conclusions are on the following two topics: 

1. Luminous AGN fraction 

We find a luminous AGN fraction of 22% ± 5% 



among massive (M* > 5 X lO^^Mg) galaxies at red- 
shifts 1.5 < z < 2.5, using their X-ray properties 
extracted from the 4 Ms Chandra Deep Field South 
observations. Among the X-ray detected galaxies, 
53% ± 13% harbor high-luminosity AGNs, while 
stacking the galaxies not detected in X-ray, leads 
to mean detections consistent with low-luminosity 
AGNs or pure star formation processes. The lu- 
minous AGN fraction among quiescent and star- 
forming galaxies is similar (19%±9% and 23%±5%, 
respectively) and does not depend on galaxy Af*. 
We confirmed that extended X-ray emission from 
a hot gaseous halo is not a viable explanation for 
the observed X-ray emission of the X-ray detected 
galaxies. 

2. Limits on total AGN fraction 

We convert the rest-frame hard band X-ray lumi- 
nosity into an upper limit on the star formation 
rate, SFR2-iokcV and compare to that derived from 
the rest-frame IR-I-UV emission, SFRir+uv- All 
luminous AGNs show an excess in SFR2-iokoV as 
expected, and so does a large fraction (~ 90%) of 
the remaining detected galaxies. While the star- 
forming galaxies not detected in X-ray have a mean 
X-ray inferred SFR of 71 ± 51M0yr~^, consis- 
tent with their SFRir+uv, the stack of quiescent 
galaxies shows an excess in SFR2-iokeV of a fac- 
tor > 10 above the upper limit on SFRir_)^uv- For 
these galaxies, we find that a minimum fraction of 
^ 60% must contain low-luminosity (-/jo.s-skeV ~ 
lO^^ergs-i) AGNs if the SFR estimates from X- 
ray are to be explained, and that low-luminosity 
AGNs might be present in all of them. On the 
other hand, for the star-forming stack, we derive a 
low- luminosity AGN fraction of 0-40%. 
Gathering all low- and high-luminosity AGNs, we 
derive a lower limit to the total AGN fraction of 
48%, with a tentative upper limit of 72%. 

Our study is the first to present observational evidence 
that, at z 2, the majority of quiescent galaxies host a 
low- to a high-luminosity AGN, while the AGN fraction 
is significantly lower in star-forming galaxies. These find- 
ings are consistent with an evolutionary scenario in which 
low-luminosity AGN quench star formation via the en- 
ergetic output from SMBH accretion, which, if believed 
to continue in an episodic fashion as often invoked by 
models, would need to have "dormant" phases at least 
as long as "active" phases. 

We find that the high-luminosity AGNs are likely re- 
lated to non-secular processes such as major mergers. In 
the future, examining the X-ray properties of galaxies 
in a larger sample, cross-correlated with signs of major 
mergers, may shed further light on the co-evolution of 
AGNs and host galaxy. 
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